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ABSTRACT 

Petrucci et al. (2002; hereafter P02) have reported extraordinary behaviour of the iron 
Ka line in the type 1 AGN Mrk 841. At the XMM-Newton/^VIC resolution, a narrow 
line was observed in a short observation, which then apparently disappeared around 
0.5 days later when another observation was performed. The limits are such that the 
line cannot both be narrow, and variable at this level, as the light-travel time places 
the material at a radius where broadening by the black hole's gravitational field is 
inevitable. Adding an additional observation taken before the P02 data, we present a 
different interpretation of the apparent variability. The data support the hypothesis 
that the line in fact varies in width, not only in flux. Rather than showing no line, 
and despite the strong Compton reflection component, the flnal observation exhibits 
broad emission at ^ 6.2 — 6.4 keV, with a width a ^ 0.8 keV, typical of a relativistic 
accretion disc. We suggest that the apparently narrow line in the early observations 
arises from local illumination by a flare inducing an hotspot in the inner disc, which 
then becomes progressively broadened as the disc rotates. 
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The fluorescence iron Kq line is a well-established feature 
in the hard X-ray spectra of Active Galactic Nuclei (AGN) 
(Nandra & Pounds 1994). The properties of this line can 
provide valuable information on the regions surrounding the 
AGN. For example, narrow Fe Ka emission can be associ- 
ated with optically thick material at large distance from 
the nuclear region (Ghisellini, Haardt & Matt 1994; Kro- 
lik, Madau & Zycki 1994) . Such narrow lines are commonly 
observed in high resolution spectra taken with the Chan- 
dra HETG (Yaqoob et al. 2001; Yaqoob & Padmanabhan 
2004). Broad lines were often observed in the ASCA spectra 
of Seyfert galaxies (Tanaka et al 1995; Nandra et al. 1997), 
and they were interpreted as the result of Doppler broad- 
ening and gravitational redshift in the proximity of a su- 
permassive black hole (Fabian et al. 1989). The existence of 
strong narrow components revealed by Chandra has caused 
the modelling of the broad emission seen by ASCA to be- 
come ambiguous in some cases, where complex absorption 
can be invoked to explain some - and in some objects all - of 
the broad residuals (e.g. Pounds & Reeves 2002, Reeves et 
al. 2004). On the other hand, the best example of a broad 
accretion disc line, MCG— 6— 30— 15, seems to survive all 
attempts at modelling the spectrum with a complex contin- 
uum (Fabian et al. 1995, 2002; Vaughan & Fabian 2004). 
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Nevertheless, given the inherent ambiguity in the spectral 
modelling, it is clearly important to seek other clues and 
indications of the origin of the emission near the iron K- 
complex. Line variability is the most obvious complement 
to the spectra. 

Mrk 841 is a type 1 AGN at z = 0.0365 (Veron-Cetty, 
2001). Several previous observations performed with differ- 
ent X-ray observatories {EXOSAT, Ginga and ROSAT) have 
reported on its spectral variability (George et al, 1993, Nan- 
dra et al 1995). Evidence for a broad iron emission line has 
been reported in one of two ASCA spectra reported by Nan- 
dra et al. (1997), and Weaver, Gelbord & Yaqoob (2001) 
have reported marginal evidence for variability of the line. 
The presence of a Compton reflection component has been 
studied and confirmed by BeppoSAX (Bianchi et al, 2001, 
2004) and XMM-Newton observations by P02. The most in- 
teresting result from the P02 observation, however, was the 
detection of apparent rapid (~0.5 days) variability of an Fe 
Kq line which was also rather narrow. P02 showed that the 
small line width observed by XMM-Newton - if due to Ke- 
plerian motions implying a distant origin for the line - was 
inconsistent with the rapid variability for any reasonable 
black hole mass (Woo & Urry 2002). This emphasises the 
ambiguity inherent in static spectral modeling as described 
above. The variability constraint in Mrk 841 is stronger than 
the spectral constraint, as the former relies only on causal- 
ity and the Eddington limit, while the latter depends on the 
uncertain velocity field of the line emitting material. 
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In this Letter an additional XMM-Newton exposure has 
been added to the one pubUshed in P02, giving better statis- 
tics and longer temporal coverage. This permits an alternar- 
tive intepretation of the P02 data, which may reconcile the 
rapid line variability with its small width. 



2 OBSERVATIONS AND DATA REDUCTION 

Mrk 841 was observed by XMM-Newton between the 13th 
and 14th of January 2001, in 3 different pointings, referred 
to as Obs I, II, III in chronological order (observations II and 
III are those published in P02) . The respective durations are 
9.3, 11.9 and 14.4 ks and I and II are contiguous, whereas 
there is a ~ 43 ks gap between II and III. The data reported 
here come from the EPIC instrument: both of the pn and 
MOS cameras were in small window mode, with only the 
M0S2 data available for analysis as the MOSl was in FAST 
UNCOMP mode. Observations II and III were performed 
with the thin filter while a medium filter was used in Obs. 
I. 

The data were reduced using the XMM-Newton Science 
Analysis Software, version 5.4.1. No background flares were 
found in the light curve, and having taken into account the 
dead time, the filtered event files contained 5.9,7.6 and 9.2 
ks of good exposure in the pn detector, for Obs I, II and 
III, respectively. Only the pn analysis is described in the 
following. The data from MOS2 have bccu analysed: they 
confirm the following results but they are not included as 
the constraints are less strong than those given by the pn. 
Source events have been collected in all 3 data sets from 
a circular region with a 40" radius centred on the source, 
while background events have been collected from a larger 
region of 56" radius, chosen in a nearby source-free portion 
of the detector. Single and double events are selected in the 
spectra yielding a number of source counts between 2-10 
keV of 9374 counts (Obs I), 13495 (Obs. II) and 17510 (Obs 
III). Pulse-height spectra were binned in order to have at 
least 20 counts/bin which allows the use of minimization 
for spectral fitting. 



3 SPECTRAL ANALYSIS 

Throughout the analysis errors are quoted at Ax^= 4.61 
(90 per cent of confidence for two interesting parameters). 
We began by fitting each spectrum individually in the 2-10 
keV rest frame range with a model consisting of a power law 

and the Galactic column density (N_ff=2.33 x lO'^" cm~^). 
The power law slopes derived for each observation were 
consistent with each other: 1.94l^-^| (I), 1.951^-24 (II) and 

i.98iro^ (in). 

Fig. 1 shows the 2-10 keV spectra plotted as the ratio of 

the observed pn data to a simple power law model, fitted to 
the spectra ignoring the 5-7 keV channels. Residuals at ~ 
6.4 keV rest-frame (dashed line) reveal the presence of an 
Fe Ka line. 

Fitting a narrow line to these individual spectra con- 
firms the result of P02, i.e. that the apparent flux drops 
in Obs. Ill If the line width is left free to vary, however, 
a different scenario is implied. As is evident from Fig. 1, 
the line in Obs. Ill does in fact appear to be present, but 




Figure 1. Data to model ratios of the three data sets (observer's 
frame), fitted with a power law with r= ~ 1.9, excluding the 5-7 
keV energy channels: the dashed line marks the rest-frame energy 
of the neutral Fe Ka line 6.4 keV. The line profile is narrow in 
the top and middle panels, while it is considerably broadened in 
the bottom one. 



is much broader than in the other two datasets. To better 
constrain this phenomenon, and since the photon indeces de- 
rived by fitting the spectra individually are consistent with 
each other, we assume the broadband spectral shape to be 
constant. Therefore, in the following we perform the spectral 
analysis for all the three observations simultaneously, allow- 
ing the cross-normalizations (and often the line parameters, 
as described) to be free. 

Adding a gaussian line with the width free to vary im- 
proves the fit with Ax^~ 74, with respect to a simple power 
law model. The rest-frame line energies are 6.29^g'4'j keV, 
6.37t.o ll keV and 6.22tg;f| keV for I, II and III respectively, 
but the line widths show some difference, with best fitting 
values 0.25, 0.13 and 0.88 keV for the three observations. 

The parameters in Obs. I and II are consistent with 
each other and these two pointings are nearly consecutive 
in time, followed by a long gap (see previous section). For 
this reason therefore, the line parameters have been tied to- 
gether combining (Obs. I + II), in order to compare with 
the broad line features seen in Obs. III. 
The best fitting parameters are summarised in Table 1. In 
Fig. 2 we plot the confidence contours of the line widths 
versus the line fluxes: it is clear from this picture that the 
primary variations is not in flux (which is consistent with a 
constant within the error bars) but the line profile - param- 
eterized here by a broadening of the gaussian width. 

It is this interesting and unusual possibility that we 
investigate further below. We have therefore continued to 
tie the parameters of (I+II) together in the fits below, and 
where appropriate have also tied parameters between Obs. 
(I-I-II) and Obs. III. In the following models the line pa- 
rameters are kept tied so to be the same in (I-I-II); also the 
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Figure 2. Confidence contours at 68.3 and 90 per cent level of 
confidence for the line widths and fiuxes when data are fitted 
with table 1 model. The left contour refers to (I+II) and the 
right contour refers to III alone. 



energies have been tied to have the same value in all the 
observations. 

The large line broadening in observation III {a ^ 0.88 
keV) can be the effect of reflection in the inner accretion 
disc which, if optically thick, will also introduce a Compton 
reflection component into the spectrum. Wc have modelled 
this using the reflection model of Magdziarz & Zdziarski 
(1995; called "PEXRAV" in the XSPEC spectral fitting 
package), with the disc inclination parameter i set to be 
30°. Wc fit this together with a neutral Fe K line, and ac- 
count for the rclativistic effects in the vicinity of the cen- 
tral black hole using the blurring code described by Fabian 
et al. (2002), in Kerr metric (Laor 1991), which we refer 
to as KDBLUR. In this way the reflection continuum gets 
smoothed and smeared and the lino broadened according to 
an emissivity law as a function of radius of the form j oc 
r~^, where j is the line emissivity and r is the radius. The 
inner and outer radii have been fixed to nn = l-24rg and 
Tout = lOOrg, respectively. The model is shown in the left 
panel of Fig. 3 and best fitting parameters arc quoted in 
Table 2. The derived emissivity index in (I-I-II) is q= — 10, 
the lowest value permitted by the model, indicating that 
the line emission in (I+II) is concentrated at large radii. On 
the other hand the line emissivity in III has q=2.1, which is 
much more typical (Nandra ct al. 1997), and indicates that 
the line is concentrated in the inner regions. The upper limit 
on the inclination angle obtained in KDBLUR implies the 
disc is consistent with being viewed at a low inclination (< 
25°). The confidence contours in the right panel of Fig. 3 
show that the refiection component is rather strong, since 
the value of the solid angle of the reflecting material seen 
by the X-ray source, R=n/27r ~ 2, is high, R = 2.llj j^. 
This is in good agreement with the data in P02, despite the 
fact that they did not incorporate the relativistic blurring, 
and that their constraints were improved by the inclusion of 
BeppoSAX high energy data (see also in Bianchi et al 2001). 
Despite the large reflection, the line equivalent widths are 
only ^ 120 cV, possibly implying a small Fo undcrabun- 
dance for the observed values of F and R in the reflection 
model of George & Fabian (1991). 

A gross change in a power-law line emissivity law seems 



an artificial and unlikely explanation for the change in line 
profile. Arguably a more plausible alternative is that the 
line arises from illumination by discrete flares acting over a 
limited radial range. A large variation in the line width could 
bo explained if the flaring regions were concentrated at large 
radii in (I-I-II) , but occurred much closer to the black hole in 
III. To test this, we set the outer radius in KDBLUR to be a 
factor of 1.1 larger than the inner radius, defining a ring on 
the disc where the line originates; here the emissivity index 
has been fixed to q= 3 and the inner radius has been left free 
to vary. With these conditions, there is a small improvement 
of Ax^=5 in the fit (see Table 3). The reflection fraction is 
slightly diminished, R ~ 1.6, better reconciling with a solar 
Fe abundance for the observed EW (between ~ 110-130 cV). 
With such a fit the emission regions must be disjointed. The 
90 per cent lower limits on the inner radii for (I-I-II) is 100 
Tg - i.e. the emission comes from outside this radius. On the 
other hand, the broad line in Obs. Ill is constrained in a very 
inner ring, at less than 20 Tg. As it is discussed in detail later, 
however, it is also possible that the line emissivity is not well 
modeled by a power law, or even circularly symmetric, so the 
results of these fits could be quite misleading. 



4 DISCUSSION 

Three XMM-Newton spectra of Mrk 841 have been analysed. 

The main results arc: 

• The 2-10 keV continuum is well fitted by a power law 
of r=2.05lS;J?. 

• The presence of a fiuorescent Fe Ka line consistent with 
neutral reflection is confirmed in all the observations. 

• Previous analysis of roughly the same data claimed ev- 
idence for variability of a factor of ~ 1.6 in the narrow line 
flux and of ~ 10-20 per cent in the continuum (P02). The 
present results indicate instead a much lower variation in the 
line flux, referring to the best fltting values, but they are still 
consistent with variability if the error bars arc taken into ac- 
count (see table 1). Anyway, a clear broadening of the line 
in observation III is detected with 99 per cent confidence. 

• The strong reflection component reported both by P02 
and Bianchi et al. (2001), is confirmed, with a value of the 
solid angle of the rcfiocting material seen by the X-ray source 
R=n/27r ~ 2. Accretion disc modelling of the broad line 
with low inclination, < 26°, and line equivalent of ~ 120 eV, 
would suggest somewhat less than solar abundance given the 
strong reflection component. 

• If modelled by a ring of emission in the accretion disc, 
the line in the two different time intervals is emitted in 
two very different regions of the putative accretion disc, one 
placed within ~ 20 r^, from the spinning black hole and one 
at > 100 Xg. 

If our interpretation of the variability - that the line 
changes in profile, rather than flux - is correct, it must 
also be reconciled with a reasonable physical interpretation 
and appropriate time scales in the disc. The key observa- 
tional facts arc that, during the first 21 ks, the lino width 
is 0.15lQ;Qg keV (I-I-II), then the detectors are switched off 
for ~ 43 ks before being operated again for the last ~ 15 
ks (III) . In this last exposure the line is broadened by Acr~ 
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Table 1. Best fitting parameters for (power law + zgauss) model with the first two observations tied together. 
The energies are quoted in the source rest-frame 

POWER LAW+ZGAUSS 



Obs. 
(Chr. order) 


r 


E 
(keV) 


(keV) 


EW 
(eV) 


Flux (line) 
(10-^ ph cm-2 B-1) 




I + II 


-0.03 


^ ..,. + 0.08 
6-*-0.06 


"■-^"'-0.08 


174+64 


, 77+0.83 
''■"-0.72 


1369/1423 


III 


idem 


idem 


n ss+0.98 

0-88-0.54 


007+246 
^■^'-162 








2 5 10 1.9 2 2.1 2.2 

Energy Pholndex 



Figure 3. Left panel: table 2 best fitting reflection model. The blurring code acts only on the line and the reflection continuum, yielding 
a very broadened line in Obs III. Right: 68.3, 90 and 99 per cent confidence contours of the photon index versus the reflection fraction 
parameter R in the same fit; the value of ~ 2.1 indicates a strong component, still consistent with the previous ones observed by 
XMM-Newton and BeppoSAX 

Table 2. Best fitting parameters for the refiection model adopted; the radii in the KDBLUR blurring routine 
are fixed to rj„=1.24 rg and rout=100rg. The line cmissivity scales with radius according to j oc r~^. 

POWER LAW+ (PEXR.W + C5AT:)*KDBLIIR 



Obs. 
(Chr. order) 


r 


R 


i 

degrees 


q 


EW 
(eV) 


X^/dof 


I + II 




^■i"-0.82 


< 26 


-10 


113+60 


1368/1423 


III 


idem 


idem 


idem 


2.08 


70+123 
'"-70 





Table 3. Best fitting parameters of the reflection model with Tout= 1-1 Xrj„ and the emissivity index q=3. 
The is by far the best obtained. The energies are in the source rest-frame. 

PO + (PEXRAV + GAU) *KDBLUR 



Obs. 


r 


R 


i 


E^ne 


EW 






(Chr. order) 






degrees 


(keV) 


(eV) 







I + II 2.05+° J^ 1.62+° |J < 26 6.41+0 56 114j:Jg > 100 1363/1420 

in i.Iciii id. 'in id. 'in i.Iciij llir, ' ■ 20 
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Figure 4. A schematic map of ttie inner accretion disc showing 
the allowed region were a narrow line with parameters consistent 
with spectrum (I+II) could be emitted. The different shades rep- 
resents regions of increasing confidence. The darkest shade corre- 
sponds to the centroid energy 6.36^g'Q4 keV, with errors within 
Ax^=l, (68 per cent of confidence for one interesting parame- 
ter). The lighter colours correspond respectively to errors of 90, 
99, 99.73 per cent. We assume the Schwarzschild case so no line 
can be generated in the innermost region within 6 Rg (lightest 
shade). 



0.73 keV. Is this consistent with any reasonable time scale 
in the system? 

Time-scales quoted in the following are in units of Mg = 
IO^Mq consistent with the black hole mass estimated by 
Woo & Urry (2002). The black hole mass was estimated 
by these authors using the Broad Line Region size derived 
from the monochromatic luminosity at 5100A in ergs s~^, 
RsLfl oc 1/5100, which combined with the BLR clouds ve- 
locity estimated by the mean Full Width Half Maximum of 
the optical spectral lines in km s~^, leads to Mbh= 4.817 

X FWHM^ Mo. We note, however, that this 

L 10 ergss ^ J ^ ' ' 

mass estimate is subject to considerable uncertainty. 

The simplest interpretation consistent with our mod- 
elling of the line is that at the time of observations (I-I-II), 
the disc was mainly illuminated at large radii, and during 
observation III, the illumination was mainly in the central 
regions, and thus broadened. A similar interpretation for 
the variable line profile in MCG--6— 30— 15 has been given 
by Iwasawa et al. (1996). 

There are, however, important differences here. Firstly, 
the difference in width is much greater, with the narrow 
line in (I+II) implying dominant emission outside 100 rg. 
It seems difficult to imagine that X-ray flares dominating 
the emission at such a large radius would then dissipate 
completely 0.5 days later. In addition, the near-constancy of 
both the line and continuum fluxes argues strongly against 
this interpretation (c.f. MCG— 6— 30— 15): there is no good 
reason why flares at such different radii would have the same 
luminosity. One could imagine a flare beginning at large radii 



and propagating inwards, retaining roughly the same power. 
However, only at the light crossing time can a flare reflected 
at ~ 100 Tg travel to ~ 10 Tg on a time-scale consistent with 
our observations. The appropriate timescale is At — 4.4 x 
10'' Mg s having assumed that the matter travels at the 
speed of light. This possibility is also unlikely: the matter in 
the disc certainly is accreted at a velocity much lower than 
c which makes the At longer than estimated above. 

Perhaps the most plausible interpretation of the data 
consists in a localized hotspot co-rotating with the disc. As- 
suming keplerian rotation, tkepi= 2.85 x 10^ Mg (r/lOra)^/^ 
s (Treves, Maraschi & Abramowicz 1988) which leads to 
tfeep!~ 1 X 10^ Mg s for a complete orbit at 10 rg, tkepi^ 
0.46 X 10^ Mg s at 6 rg and ~ 0.043 x lO'^ Mg s at the last 
stable orbit r= 1.24 rg. Thus, the rotation time is consistent 
with the variability timescale as long as the action occurs at 
a small radius. 

An isolated flare on the disc would illuminate a nar- 
row range of azimuthal angle (and hence Doppler factor), 
resulting in a narrow line. As demonstrated in Dovciak et 
al. (2004), the blue horn originated by an internal spot (6 rg) 
in a disc observed at 30°, is very close to 6.4 keV and it can 
easily be misinterpreted as a narrow line originating further 
away. These authors calculate and display line profiles inte- 
grated over 1/12 of the orbital period at the corresponding 
radius. In a disc viewed at 30°, a line produced at r=10 and 
20 rg can result in a narrow peak at almost 6.4 keV. The line 
in (I+II) is consistent with this scenario, being integrated 
over 1/10 of the estimated orbital time at ~25 ° and 20 

Over the Keplerian timescale, as the disc rotates, the 
line could be broadened out as seen in observation HI. The 
physical mechanism for the broadening mechanism and is re- 
lationship to the rotation is unclear, but it must be such that 
the hotspot is well confined during the first part of the disk 
orbit and subsequently, the emissivity is distributed over a 
larger region on the disc. It is also possible that a shift of 
the spot inwards is contributing to produce a broader line, 
even though it is unfortunately very difficult to place tight 
constraints on the size or the position of the flaring region. 
The model proposed by Dovciak et al. predicts that both 
line energy and flux should present some amount of varia- 
tion, over the orbital timescale; such variations cannot be 
tracked and measured precisely from the available spectra, 
but by means of a qualitative comparison with fig. 2 in Dov- 
ciak et al., we can say that the AE and AF inferred from 
our observation, are consistent with what is theoretically ex- 
pected. This is perhaps unsurprising given the constraints: 
for example, when left free, the line energy in Obs HI spans 
a wide interval (6.25lo'7i keV). It is therefore easily possi- 
ble to find a region of parameter space consistent with the 
present observation as long as the spot is located in the inner 
disc during the last part of the orbit, as defined by Dovciak 
in fig. 2 (i.e. Time = 0.7-1.0). 

The phenomena discussed here are probably strongly re- 
lated to the transient, narrow and shifted emission lines now 
reported for a number of XMM-Newton and Chandra spec- 
tra of AGN (Turner et al 2002, Guainazzi 2003, Yaqoob et 
al. 2003, Turner et al. 2004, Porquet et al. 2004 ). The main 
difference in the case of Mrk 841 is that we see an apparent 
evolution of the line profile from narrow to broad. Although 
we cannot strongly rule out the possibility with the present 
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data that the broad component is actually present in the 
Obs. (I+II) spectrum, (the 90 per cent upper limit on the 
broad line in (I+II) is 162 eV), our interpretation of a pro- 
file variations seems the clearest and most straightforward 
interpretation of the data. 

There are though some effects to be taken in account 
which require accurate tuning. When the line fluorescence 
takes place so close to the black hole, as it is case, the Fe Ka 
photons are subject to a substantial gravitational redshift 
according to E' /Ea= (l-2rg/r)^/^ (Fabian et al. 1989). At, 
e.g., 20 Tg the line centroid would be shifted from the rest 
energy 6.4 keV to 6.07 keV. The Doppler shift must also be 
taken into account, and in particular if the flare originates 
in the approaching part of the disc it will be blueshifted, 
counterbalancing the gravitational redshift and, potentially, 
resulting in a narrow line close to 6.4 keV, as observed in 
Obs. (I+II). 

To illustrate these effects, we have left the line energies 
untied in (I+II) and III and calculated the appropriate grav- 
itational and Doppler shifts in a Keplerian accretion disc for 
(I+II). Then, we constructed a probability map of the ex- 
pected line energy for each radius and azimuthal angle on 
the disc. The map is displayed in Fig. 4: the maximum radius 
is 20 Tg as we determine by fitting spectrum III and we con- 
sider the Schwarzschild case, so that no line originates within 
6 Tg (lightest shade). The line energy is computed following 
the approximations in Fabian et al. (1989) and assuming the 
best fitting values in Table 3 for the disc inclination and the 
upper limit on the emission radius. The confidence regions 
refer to the errors on line energy quoted at Ax'^= 1, 2.71, 
6.63 and 9 for one interesting parameter, corresponding to 
68, 90, 99 and 99.73 per cent of confidence respectively. The 
energy ranges are illustrated in colour scale, starting from 
the 68 per cent region in the darkest shade (see caption for 
details) . 

The map demonstrates that the narrow line in Obs. 
(I+II) is consistent with being emitted within the same ra- 
dius found for the line in Obs. Ill, although from a restricted 
area on the inner disc. It is computed in Schwarzschild met- 
ric excluding the region within 6rg. However as the allowed 
region is outside this radius it should be fully consistent with 
the case of spinning black hole. 

A final caveat is to remember that the masses estimated 
from optical luminosity suffer considerable uncertainty, as 
Woo & Urry warned themeselves. The constraints on the 
emission region, and indeed on the line variability overall, 
would be loosened considerably if this were taken into ac- 
count. A prediction of our model is that the line seen in 
Mrk 841 during its "narrow" phase is not truly narrow, and 
should have a finite and detectable width - albeit not at the 
XMM/EPIC resolution. New observations at high resolution 
with, e.g. Chandra/RETG or ASTROE-H wiU be required 
to measure the width of the line more precisely, and distin- 
guish definitively between the various options. 
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